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Abstract: We examine the effects of QCD initial-state radiation on a class of quantities, 
designed to probe the mass scale of new physics at hadron colliders, which involve longitu- 
dinal as well as transverse final-state momenta. In particular, we derive universal functions 
that relate the invariant mass and energy distribution of the visible part of the final state 
to that of the underlying hard subprocess. Knowledge of this relationship may assist in 
checking hypotheses about new processes, by providing additional information about their 
scales. We compare our results with those of Monte Carlo studies and find good general 
agreement. 
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1. Introduction 

Searching for new piiysics at liadron colliders is difficult, partly because of the complexity 
of the expected new signals (typically, multiple jets and/or leptons plus missing energy) 
but also because processes at high mass scales are accompanied by copious initial-state 
QCD radiation (ISR). As this tends to be emitted close to the incoming beam direc- 
tions, and the longitudinal momentum of the hard subprocess is anyway unknown, search 
variables [1-41] have generally been constructed either from the transverse components 
of observed final-state momenta or else by assuming that some subset of these momenta 
can be unambiguously ascribed to the hard subprocess. Variables that make use of all 
observed momenta without hypothesizing any particular structure of the final state are 
termed global inclusive variables in ref. [41]. Examples of transverse global quantities of 
this type are the observed transverse energy Et, the missing transverse energy and 
their sum Ht = Et+ ^t- The distributions of such quantities can provide information on 
the energy scales of new processes such as super symmetric particle production [1,2,6]. 

Although the longitudinal components of final-state momenta are strongly influenced 
by ISR, they do contain information about the underlying hard subprocess. Indeed, the 
amount of ISR emitted is determined by the energy scale of the subprocess. It is therefore 
of interest to quantify the effects of ISR on global observables that involve longitudinal 
momentum components. The aim of the present paper is to take the first steps in that 
direction. 

In ref. [41] various global search variables were investigated, including those that make 
use of longitudinal as well as transverse momentum components. The quantities studied 
included the total energy E visible in the detector and the visible invariant mass M, 




(1.1) 
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where Pz is the visible longitudinal momentum. In addition a new variable was introduced, 
defined as 

where the parameter Mjnv is a variable estimating the sum of masses of all invisible particles 
in the event: 



Mi^^ = Y,mi- (1.3) 



i=l 

It was argued that the peak in the distribution of sj^-^^ is a good indicator of the mass scale 
of new physics processes involving heavy particle production. 

It will also be useful to define the rapidity of the visible system. 

The present paper examines the effects of ISR on global inclusive search variables, first 
in an approximate fixed-order treatment taking into account collinear-enhanced terms, and 
then in an all-orders resummation of such terms. We quantify the way the distributions 
of quantities that involve longitudinal momenta depend on the scale of the underlying 
hard subprocess and on the properties of the detector, in particular the maximum visible 
pseudorapidity ??max- With the insight thus gained, it may be possible to correct for this 
dependence and thereby extract information on the hard subprocess from such quantities. 

2. Fixed-order analysis 

The Monte Carlo results presented in ref. [41] show that the second term on the right-hand 
side of eq. (|l.2|) is not strongly affected by ISR. The first term is intended to add extra 
longitudinal information about the hard subprocess, allowing a more reliable determination 
of its mass scale. The extra longitudinal information enters through the visible mass M. 
We therefore concentrate mainly on this quantity. 

2.1 Born approximation 

Since the effect of invisible final-state particles is not the central issue here, let us suppose 
first that all the final-state particles from the hard subprocess are detected.^ Then in Born 
approximation, assuming that no beam remnants are detected, M yields a perfect estimate 
of the centre-of-mass energy of the hard subprocess. For incoming partons with momentum 
fractions xi.2, 

E=^VSixi+X2) , Pz = ^VSixi - X2) , (2.1) 

where x/S" is the hadron-hadron centre-of-mass energy, so that 

M^ = xiX2S, y = -ln— . (2.2) 

2 X2 



^We comment in sect. 5 on the treatment of invisible particles. 



- 2 - 



The differential cross section for parton flavours a, b is thus 

^^""^ ^ dxi dx2 faixi)fb{x2)S{M'^ - xiX2S)5 ( y - In — ) aab{xiX2S) (2.3) 



where fa,b are the relevant parton distribution functions for the incoming hadrons and aab 
is the hard subprocess cross section. Hence at Born level we find 



The parton distributions are normally given as Fi{x) = xfi{x), in terms of which we 
have 

.2 daab y\ ^ [ M 



" MiiiY = ^' Its' j ' ■ 

If the partonic cross section dab has a threshold or peak, indicating that the ah subprocess 
has a characteristic scale Q, then this is also manifest in the Born cross section ( p.5| ) at 
M Q, provided the relevant parton distributions are large enough for that subprocess to 
contribute significantly. 

From eqs. and ( |2.2D we have E = Mcoshy and therefore the visible energy 

distribution is given in Born approximation by 



^2 daab ^ j^2 daab 



(2.6) 

M=E sech Y 



dE'^dY dW^dY 
2.2 Quasi-collinear NLO correction 

To examine the sensitivity of the above results to ISR, let us first compute the NLO 
contribution due to quasi-collinear gluon emission and the associated virtual corrections. 
Consider first the emission of a gluon from parton a. If the emission angle 6 is large 
enough, say 9 > 6^, the gluon enters the detector and contributes to M. In the small-angle 
approximation we then have 

E = ^VS{xi/z + X2) , P, = ^^{xi/z-X2) , (2.7) 
where xi/z is the momentum fraction of parton a before the emission, so that 

M^ = xiX2S/z, y = -ln— . (2.8) 

2 ZX2 

The correction associated with a detected emission from parton a is then 

as f d9 dz ^ ^ fS / n r / n r / \r/]\^2 ^ o/ -'^1 xi \ 



dxidX2Pa{z)fa{xi/z)fb{x2)d{M - X1X2S/ z)6 [ Y - -In ] aab{xiX2S) 

TT Jq^ e Z \ 2 ZX2J 

(2.9) 

where Pa{z) is the unregularized a ^ ag splitting function. 

On the other hand if the gluon misses the detector [9 < 9c), E and Pz are still given 
by (2.1), so the contribution is 

— / ^ — dxidx2Pa{z)fa{xi/z)fb{x2)6{M^ - XiX2S)5 [Y -l;ln—] aab{xiX2S) . 
TT Jo 9 Z \ 2 X2j 

(2.10) 
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Finally the associated virtual correction is the term that regularizes the splitting function, 
which in this approximation is 



-?/ 



de 



dz dxi dx2Pa{z)fa{xi)fh{x2)5{M'^ - XiX2S)5 [ Y 



iln^ 

2 X2 



aabixiX2S) . 

(2.11) 



Adding everything together gives a correction 



/ _dOab_\ _ as f d0_ 



\dM'^dY ) TT J e 



■ dz dxi dX2 Pa{z) fh{x2)cFah{xiX2S) 



-fa{xi/z)6 (y-lln ^) 6{M^ - xiX2S/z)e{9 - 6,) (2.12) 

Z \ 2 ZX2 J 

+{-Ja{xi/z)e{e, -6)- fa{xi)]5 {y-\ In ^"j 5{M^ - X1X2S) 



X2J 

Setting aside for the moment the possibility of splittings other than a ag, the 
DGLAP evolution equation for fa{xi) is 



l-^fa{xi) = I dzPa{z) 



^fa{xi/z) - fa{xi) 



(2.13) 



where q represents the scale at which the parton distribution is measured. Hence the 
correction may be written as 



'b{x2)^ab{xiX2S) 



+ 



as f dz 



TT 



dq \ 2 X2 J 
1 , Xl 



X2S) 



'-Pa{z) fa{xi/z) <! M ^ - § In ^ j 5{M^ - X1X2S/Z) 



-6[y 5{M'^ - X1X2S) 1 6(0 - 



(2.14) 



Since d9/6 = dq/q, the first term represents a change of scale in the Born term. It replaces 
the reference scale in fa by the scale Q of the hard subprocess. The remaining terms give 
a correction 

M 



f duab \ _ as _ 
[dM^dYj ~ TrSJe^ 9 



d9 f dz ~ 

—Pa{z)fb 

z 



S 



e ^ ] X 



fa 



M 

7S' 



Z&abizM^) - fa 



M 



Oab{M^ 



(2.15) 



In leading-log approximation the 9 integration just gives a factor of — In 0c. In the same 
approximation, we may set — In^c = ??max, the maximum pseudorapidity seen by the de- 
tector. Note that this is a different quantity from Y , the true rapidity of the visible system, 
given by eq. ( |1.4[ ). The correction associated with parton h gives the same expression with 
a ^ h and Y ^ —Y . Thus, defining 



My- M _Y 

Xl = —;=€ , X2 = —;=e 
^/S y/S 



(2.16) 
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Figure 1: Distribution of visible mass M in tt production at LHC for ?7„iax = 5 and Y — 0: leading 
and approximate next-to-leading order. 



we have 



S 



d(Tab 

dM^dY 



faixi,Q)fb{x2,Q)aab{M^ 



+ Vn 



as f dz 



IT 



Z L 



z{Paiz) + Pbiz)}fa{xi,Q)Mx2, Q)aabizM^) 



- {Pa{z)fa{xi/z,Q)h{x2,Q)+Pb{z)fa{xi,Q)h{x2/z,Q)]aab{M^) (2.17) 



Expressing this in terms of Fi{x) = xfi{x), as in eq. ( |2.5| 
da„h 



dJVPdY 



Fa{xi,Q)Fi,{x2,Q)aab{M') 



+ r/max— / dz 
IT 



{Paiz) + Pb{z)}Fa{xi, Q)Fb{x2, Q)aab{zM^) (2.18) 
{Pa{z)Faixi/z, Q)Fbix2, Q) + Pb{z)Fa{xi,Q)Fb{x2/ z, Q)]aab{M'') 



Results for ti production at the LHC {pp at ^fS = 14 TeV) with r/max = 5 and 1" = are 
shown in fig. |^. Leading-order MSTW parton distributions [42] were used. For simplicity 
we have taken Q = M . Recall that the simplifying assumption made here is that all ti 
decay products are detected, so the M distribution vanishes below ti threshold. We see 
that there is a large negative NLO correction near threshold, followed by a broad positive 
peak. 

To understand these qualitative features, consider the case a = 6, as in gg — > tt, and 
y = 0, so that xi.2 = M/y/S = x. Then the NLO correction becomes simply 



5 



da 



dM^dY / 



2r,m.^^F{x) I dzPa{z) [F{x)a{zM^) - F{x/z)a{M^)\ (2.19) 
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The first term is positive-definite, contributes only above thresliold, and diverges at thresh- 
old. It produces the broad positive peak. The second term is negative-definite, contributes 
around threshold, and has a divergent coefficient. It provides the sharp negative peak. 

From eqs. ( |2.7[ ) and ( |2.8D , the relation E = Mcoshy still holds in this approximation 
and therefore the visible energy distribution is again given by eq. ( |2.6| ). 



3. Resummation 

By adding and subtracting the expression 

{Pa{z) + Pb{z)}Fa{xi,Q)Fb{x2, Q)Oab{M^) (3.1) 



in the integrand of eq. ( 2.18 ) and comparing with eq. ( |2.13 ), we see that the last line of 



that equation corresponds to a change of scale Q — > Qc = QcQ in the parton distributions, 
leading to 



M^^^ = Fa{xiMF,{x2M^ab{M^) (3.2) 



where to first order 



+ r/max^ J dz{Pa{z) + Pk{z)}{aab{zM^) " ^abiM^)} ■ (3.3) 

The interpretation of this result is simple: undetected ISR at angles less than 9c, corre- 
sponding to scales less than 9cQ, is absorbed into the structure of the incoming hadrons. 
To resum the effects of gluons at angles greater than 9c, consider first the real emission of 
n such gluons from parton a. In the quasi-collinear approximation these form an angular- 
orders sequence, giving rise to a contribution to T,ab of 

as\" f d9i f d92 f d9n ^_ /. n fs /.. /.. ., ,^2 



- - I dzi . . . dZnPa{zi) ■ ■ ■ Pa{Zn)^ab{zi ■ ■ ■ ZnM ) 

71" ^ Jec ^1 Jei ^2 Jd„.i Jo 
1 / as\'^ ~ ~ 2 

= -r ??max dzi... dZnPa{zi) ■ ■ ■ Pa{Zn)^ab{zi ■ ■ ■ ZnM ) , (3.4) 

n! V TT J Jo 

where again we have made the identification — In^^ = ^max- The multiple convolution of 
the momentum fractions Zi can be transformed into a product by taking moments. Defining 

CO 

dM^{MY'' aab{M^)^a-^ (3.5) 



we have 

( ?)" "^^^ (M^)""^ j dzi... dzMzi) . . . Pa{zn)a{zi ■ ■ ■ z^M^) = {TnT (3.6) 

where ^ 

Tn = — f dzz^'-'Paiz) . (3.7) 



vr 
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Therefore defining correspondingly 

fOO 

JO 

the contribution (|3.4|) to this quantity will be 



1 



-(r/max7^)"a^^ (3.9) 
n\ 

which summed over n gives 

exp (r?max7^) <7^^ (3.10) 
The corresponding virtual contributions give a Sudakov-like form factor 

and therefore the total contribution from parton a is 

exp (r?i„ax7^) (3-12) 
where 7^ is the anomalous dimension 

1% = — C dz (z^-l - 1) Pa{z) = ^ [\z z'^-'Paiz) , (3.13) 
Jo Jo 

Pa{z) being the regularized a ^ ag splitting function. Parton b gives a similar factor with 
7^ in place of 7^, so the result for the quantity (3^) is simply 



j^ab ^ gW(7^+7?,)^^fe . (3.14) 



We can see as follows that this result is qualitatively correct. The anomalous dimensions 

■'N 



are positive for small N and negative for large A^. Thus, for 6c <C 1, is enhanced relative 



to fj^ at small and suppressed at large N. Now from the moment definition ( |3.5| ) small 
N corresponds to large M and vice versa. Hence the distribution of M is suppressed at 
small M and enhanced at large M relative to the Born term, as observed in the Monte 
Carlo [41] and NLO results. 

The emission of partons other than gluons is included by introducing the anomalous 
dimension matrix T^r with elements given by 

{rN)ba = — C dzZ^'-^PUz) (3.15) 



where Pba{z) is the regularized a ^ h splitting function. Then 

= (e"— ^^),,, [e''-''^''),,, • (3.16) 
The corresponding generalization of the evolution equation ( 2.13| ) is 



q^Mx) = ^ / '^PUz)fa{x/z) (3.17) 
oq TT J z 
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Defining the moments of the parton distribution functions 

r.1 



rN 



dxx^'-'faix) 







we see that 



with solution 



d 



fNil) = {[Q/Qof''),JUqo) . 



Hence 



where 



fMQ) 



'')bafN{Qc 



Qc = OcQ = Qe' 



(3.18) 

(3.19) 

(3.20) 
(3.21) 
(3.22) 



showing that the evolution of the visible mass distribution is related to that of the parton 
distributions over the same range of scales. 

Taking into account the running of the strong coupling as{q) in the evolution equation 
( 3.17 ), eq. ( 3.21 ) becomes 



where 



with p = 6/(11Ca — 2nj) and 



T^ba 



vr 



Qs(Qc) 
asiQ) 



ba 



{'^N)ba = {^N)ba 

as 



I'dzz'^-'PUz) 
Jo 



The running of as will affect eq. ( 3.16| ) similarly, giving 



Sab ;^a,'b' r^a'a j^b'b 
N — CFn 



This implies that 



where 



^ab{M^)= I dzaa'b'{zM^)Ha>b',ab{^) 

Jo 



dzz^ ^Ha'b',ab{z) 



or, inverting the Mellin transformation, 

Ha'b',ab{z) 



2m 



dNz-^Ki'^Kf 



(3.23) 
(3.24) 

(3.25) 

(3.26) 
(3.27) 
(3.28) 

(3.29) 



c 



where the contour C is to the right of all singularities of the integrand. 
Alternatively, eq. ( 3.27 ) can be expressed as a double convolution, 

llabiM"^) = / dZidZ2Cra'b'{ziZ2M'^) Ka'a{zi) Kyb{z2) 

Jo 



(3.30) 
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where 

Kf= I dzz''-^K,,f,iz) . (3.31) 



It then follows from eq. ( 3.24 ) that Kyi,{z) obeys an evolution equation hke that of the 
parton distributions: 

Q^Ky,{z) = / ^P,,^{z')K,,{z/z') . (3.32) 

OQ TT J Z' 

Putting everything together, the visible mass distribution is related to the hard sub- 
process cross section (in the absence of invisible final-state particles) as follows: 

dMHY ^ j dzidZ2aa'b'{ziZ2M^) Ka>aizi)Fa{xi,Qc) Kh'b{z2)Fbix2,Qc) (3.33) 

where the kernel functions Ka'a{z) and Ky^iz) can be obtained by solving the evolution 
equation ( 3.32| ) with the initial condition that Kab{z) = 5ab5{^ — z) sA, Q = Qc- 

Since the partons sampled at scale Qc are always regarded as (anti-)collinear, the 
relation E = M cosh Y still holds and therefore the visible energy distribution is given by 



^2 dOgb ^ ^2 dcTab 



dE'^dV dM'^dV 



(3.34) 

M=EscchY 



in leading-logarithmic approximation to all orders. 

To verify that the integrated cross section is not affected by resummation, define 
2^1,2 = -21,22^1,2 and write eq. ( p. 33 ) as 



^2 daab 



M ^j^2^Y ^ ' '^^^'^^'^^'^'b'i^^^'^'^) ^a'a{xi/xi)fa{xi,Qc) Kb'b{x2/x2)fb{x2,Qc) ■ 



Now 



and 



Hence 



(3.35) 

^dY = ^^ (3.36) 

M2 Xi X2 

^ I ^Ka'a{xi/xi)fa{xi,Qc) = fa'{xi,Q) ■ (3.37) 

^ / dM^dY = ^ / dxidx2^a'b'ixiX2S)fa'ixi,Q)fb'ix2,Q) (3.38) 

ab a'b' 

in agreement with eq. ( |2.3D . 

Resummed results corresponding to fig. |l] are shown in fig. We see that the peak of 
the distribution has moved to much higher mass, beyond 1 TeV. This is due to multiple 
emission of ISR partons in the evolution of the initial state from the detection scale Qc 
to the hard subprocess scale Q. As the value of 7?max is reduced, the range of evolution 
becomes smaller, less ISR is emitted, and the peak moves closer to the hard subprocess 
scale, as illustrated in fig. H. However, recall that any loss of visible particles produced 
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Figure 3: Resummed distribution of visible mass M in tt production at LHC iorY — Q and lower 
values of rymax^ colour scheme as in figs.l and 2. 

in the hard subprocess (top decay products in this case) has been neglected here. Such 
losses will cause the peak to fall below the hard subprocess scale at low values of 7?max- In 
ref. [41] it was found that when r/max = 1-4 the peak lies close to threshold for the hard 
subprocesses studied there, presumably as a result of compensation between ISR and loss 
of visible decay products. 

Results for higher values of the visible rapidity Y are shown in fig. ^. The peak moves 
to lower mass as Y increases, as a consequence of the suppression of high masses by the 
rapid fall-off of the parton distributions at high x. 

4. Monte Carlo comparisons 

In this section we compare the predictions of the above analytical treatment with Monte 
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Figure 4: Resummed distribution of visible mass M in tt production at LHC for Tymax = 5: results 
at non-zero visible rapidity Y. 




Figure 5: Monte Carlo results for various inclusive observables in tt production at LHC. The yellow 
filled histogram shows the true hard subprocess invariant mass. The dashed red and dot-dashed 
blue curves show the results of eqs. (3.33) and (3.34) , respectively, integrated with respect to Y. 



Carlo results from HERWIG [43,44]. Figure ^ shows the Monte Carlo results for various 
global inclusive observables in tt production at the LHC when r/max = 5. For the detector 
simulation we used GETJET [45] with calorimeter cell size (Ary, Ac/)) = (0.1,0.1). To 
facilitate the study of ISR effects, the simulated underlying event was turned off. 

As expected, the distributions of variables that involve longitudinal momenta are 
broadened and shifted relative to those of purely transverse quantities. The distributions 
of the visible mass M and the new quantity (1.2) are similar, while the visible energy E 
has a broader distribution, as it includes the visible momentum. 

Comparing the visible mass and energy distributions with the analytical calculations 
(the dashed red and dot-dashed blue curves, respectively), we find good overall agreement. 
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Figure 6: Position of peaks in the distributions of various inclusive observables in ti production at 
LHC, as a function of the maximum visible pseudorapidity rymax- Solid curves: Monte Carlo results. 
Dashed red and dot-dashed blue curves: M and E peaks from eqs. ( 3.33 ) and ( 3.34 ), respectively. 
Dotted line: tt invariant mass peak. 



considering the simplifications made in the latter, viz. the quasi-collinear approximation, 
no loss of top decay products and no hadronization. In addition, turning off the underlying 
event in the Monte Carlo does not entirely eliminate contributions from spectator parton 
fragments at high rapidity. 

The motion of the peaks in these distributions as the maximum visible pseudorapidity 
Vmax is varied is shown in fig. |^. The distributions of the transverse quantities remain 
remarkably stable, while the peak positions of the other variables rise faster than linearly 
with increasing ??max- The analytical prediction for the peak position of the visible mass 
is in good agreement with the Monte Carlo results at intermediate values of ?ymax 

but 

somewhat higher at low values and lower at high values. This is consistent with the loss 
of top decay products at low ??max and the contribution of spectator fragments at high 
^max in the Monte Carlo. The peak position of the visible energy distribution is somewhat 
overestimated at all but the highest values of ?7max- However, as may be seen from fig. |5|, the 
distribution of this quantity is quite flat near the peak, and therefore a small discrepancy 
between the analytical and Monte Carlo results can give rise to a larger difference in the 
associated peak positions. 

More precise Monte Carlo studies of global inclusive observables like those considered 
here could be performed using methods that match parton showers to matrix elements 
for extra jet production; see for example [46-48] and references therein. We checked that 
matching to NLO (up to one extra jet) using MC@NLO [49,50] did not lead to significantly 
different results. 
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5. Conclusions 



We have derived a simple analytical formula, eq. ( 3.33| ), for the effects of QCD initial state 
radiation on the invariant mass M and rapidity Y of the visible part of the final state 
of a hadronic collision, as a function of the pseudorapidity coverage ?7max of an idealized 
detector. Given the mass distribution of a hard subprocess involving incoming partons a' 
and 6', a simple convolution with the universal functions Ha'h',ah in eq. ( p. 29]) provides the 
quantities E^fe that, when multiplied by the relevant parton distributions, sampled at an 
^/max-dependent scale, yield the visible mass and energy distributions. 

The derivation involved a number of simplifying assumptions and approximations, but 
was in satisfactory agreement with Monte Carlo results that should in principle be more 
realistic. The difference between the visible mass M and the new variable ( |1.2D is small 
when Minv = 0, essentially because of the relative smallness of the missing transverse 
energy manifest in fig. ^. 

Many of the approximations made here could be improved if a more precise analysis 
is needed. The emission of invisible particles in the hard subprocess could be taken into 
account by replacing the mass distribution of that process by a visible mass distribution, 
with the invisible component already integrated out. The universal functions Ha'b' ,ab could 
be computed to next-to-leading order, as they involve the same quantities that drive the 
evolution of the parton distribution functions. Together with an NLO calculation of the 
hard subprocess mass distribution, this would provide a complete NLO description of the 
visible mass and energy distributions. 

Whether such refinements are worthwhile depends on the extent to which observables 
involving longitudinal momenta are found useful in the exploration of physics beyond the 
Standard Model. The Monte Carlo results presented in the previous section confirm that 
transverse quantities are much less sensitive to the effects of ISR. However, at the very least 
it will be useful to check the consistency of hypotheses about new subprocesses with the 
distributions discussed here, which do contain independent information about their scales. 
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